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The effect of a gate voltage (V g ) on the spin-splitting of an electronic level in a quantum dot 
(QD) attached to ferromagnetic leads is studied in the Kondo regime using a generalized numerical 
renormalization group technique. We find that the V^-dependence of the QD level spin-splitting 
strongly depends on the shape of the density of states (DOS). For one class of DOS shapes there 
is nearly no V^-dependence, for another, V g can be used to control the magnitude and sign of the 
spin-splitting, which can be interpreted as a local exchange magnetic field. We find that the spin- 
splitting acquires a new type of logarithmic divergence. We give an analytical explanation for our 
numerical results and explain how they arise due to spin-dependent charge fluctuations. 
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The manipulation of magnetization and spin is one 
of the fundamental processes in magneto-electronics and 
spintronics, providing the possibility of writing informa- 
tion in a magnetic memory |l[ , and also because of the 
possibility of classical or quantum computation using 
spin. In most situations this is realized by means of an 
externally applied, nonlocal magnetic field which is usu- 
ally difficult to insert into an integrated circuit. Recently, 
it was proposed to control the magnetic properties, such 
as the Curie temperature of ferromagnetic semiconduc- 
tors, by means of an electric field: In gated structures 
0, due to the modification of carrier-density- mediated 
magnetic interactions, such properties can be modified 
by a gate voltage. In this Letter we propose to control 
the amplitude and sign of the spin-splitting of a quantum 
dot (QD) induced by the presence of ferromagnetic leads, 
only by using a gate voltage without further assistance of 
a magnetic field. To illustrate this effect we investigate 
the Kondo effect and its spin-splitting as a very sensitive 
probe of the spin state of the dot and the effective local 
magnetic field in the QD generated by exchange interac- 
tion with the ferromagnetic leads. 

Recently, the possibility of the Kondo effect in a QD at- 
tached to ferromagnetic electrodes was widely discussed 
d II H E3 II la, and it was shown, that the Kondo 
resonance is split and suppressed in the presence of fer- 
romagnetic leads [3, 13 . It was shown that this splitting 
can be compensated by an appropriately tuned external 
magnetic field, and the Kondo effect is thereby restored 
0, Q • In all previous studies of QDs attached to ferro- 
magnetic leads 0, 0, 0, |y, 0, 01 IaJ an idealized, flat, spin- 
independent DOS with spin-dependent tunneling ampli- 
tudes was considered. However, since the spin-splitting 



arises from renormalization effects i.e. is a many-body ef- 
fect, it depends on the full DOS-structure of the involved 
material, and not only on its value at the Fermi sur- 
face. In realistic ferromagnetic systems, the DOS shape 
is strongly asymmetric due to the Stoner splitting and 
the different hybridization between the electronic bands 

Q. 

In this Letter we demonstrate that the gate voltage de- 
pendence of the spin-splitting of a QD level, resulting in a 
splitting and suppression of the Kondo resonance, is de- 
termined by the DOS structure and can lead to crucially 
different behaviours. We apply the numerical renormal- 
ization group (NRG) technique extended to handle bands 
of arbitrary shape. For one class of DOS-shapes, we find 
almost no V^-dependence of the spin-splitting, while for 
another class the induced spin-splitting, which can be in- 
terpreted as the effect of a local exchange field, can be 
controlled by V g . The spin-splitting can be fully compen- 
sated and its direction can even be reversed within this 
class. We explain the physical mechanism that leads to 
this behavior, which is related to the compensation of the 
renormalization of the spin-dependent QD levels induced 
by the electron-like and hole-like quantum charge fluc- 
tuations. Moreover we find that for the QD level close 
to the Fermi surface, the amplitude of the spin-splitting 
has a logarithmic divergence, indicating the many-body 
character of this phenomenon. 

Model and method. - The Anderson model (AM) of a 
single level QD with energy eo and Coulomb interaction 
U, coupled to ferromagnetic leads, is given by 
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Here c r ka and c^- (n CT = d) a d a ) are Fermi operators for 
electrons with momentum k and spin a in the leads 
(r = L/R), and in the QD, V^ is the tunneling am- 
plitude, S z = (n-f — fii)/2, and the last term denotes the 
Zeeman energy of the dot. The energy eo is experimen- 
tally controllable by V g (eo — V g ). 

In order to discuss the gate voltage dependence of the 
QD level spin-splitting, we consider here a more real- 
istic, both energy and spin dependent band structure 
[p rT (cj) ^ Prl(w)], violating p-h symmetry p r<T (u) ^ 
p ra (—u), which leads to an energy dependent hybridiza- 
tion function T r(T {uo) = 7r^2 k 8{uo-e ka )V^ k = 7ip ra (uo)V^, 
where we take V r k = V r to be constant. We apply the 
NRG method |l(l.lllj extended to handle arbitrary DOS 
shapes and asymmetry. To this end, the standard log- 
arithmic discretization of the conduction band is per- 
formed for each spin component separately, with the 
bandwidths, D^ = D^ — Do, chosen such that the to- 
tal spectral weight is included in [—Do, Do] for all values 
of V g studied here (to avoid different systematic errors 
upon changing V g ). 

Within each interval [— uo ni — a; n +i] and [u; n +i , uj n ] 
(with uo n = DoA~ n ) of the logarithmically discretized 
conduction band (CB) the operators of the continuous 
CB are expressed in terms of a Fourier series. Even 
though we allow for a non-constant conduction electron 
DOS, it is still possible to transform the Hamiltonian 
such that the impurity couples only to the zeroth order 
component of the Fourier expansion of each interval |l2| . 
Dropping the non-constant Fourier-components of each 
interval [lQ, LUj then results in a discretized version of 
the Anderson model with the continous spectrum in each 
interval replaced by a single fermionic degree of freedom 
(independently for both spin directions). Since we allow 
for an arbitrary DOS for each spin component a (|, j) of 
the CB this mapping needs to be performed for each a 
separately. This leads to the Hamiltonian: 



H = ^ e °™° + U ^?ll + V^Oa/TT^KJr/Oa + fL d *] 
a a 

oo 

i / J [^naJnaJna + tna\J na Jn+la + Jn+lcrJncr)\ •> (*) 
an=0 

where f ncr are fermionic operators at the nth site of the 
Wilson chain, £oa = 1/2 J_ D ° T a (uj)duj J t na denotes the 
hopping matrix elements, and e a — eo — BS Z . The ab- 
sence of particle-hole symmetry leads to the appearance 
of non-zero on-site energies, e na along the chain. In this 
general case no closed expression for the matrix elements 
t na and £ ncr , both depending on the particular structure 
of the DOS via T a (u), is known, therefore they have to 
be determined recursively. This requires rather advanced 
numerical methods, due to the exponentially fast decay 
of t ncr and e ncr along the chain 
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FIG. 1: (color online). V g -dependence of the spin- splitting: 
Normalized spectral function tt ^2^ r a (0)A a (oo) as a function 
of energy cj and gate voltage eo, for the three different DOS 
shapes (depicted in insets) characterized by a different Q, 
which modifies both the spin and p-h asymmetry: (a-c) for 
magnetic filed B = 0, (d) B/U = 0.017, (e), and (f) B/U = 
0.0083. The white dashed lines are obtained using Eq. @. 
Here U = 0.12L> , ttV<? = UD/6, A = 0.15D and T = 0. 
Inset: the scheme of the parabolic DOS shape for spin | (red) 
and | (blue). 



This method allows one to calculate the level occu- 
pation Ha = (her) and the spin-resolved single-particle 
spectral density A a (u) = — -ImQ v a ((jS) , where G v a {w) de- 
notes a retarded Green's function. For symmetric cou- 
pling [Tl(t(w) = TR a (uj)] the spin-resolved conductance 
takes the form G a = tt^- J_^ duT (T (uj)A (7 (uj)( — q^~) 
where f(uo) is the Fermi function. 

Spectral function and conductance. - Here, we fo- 
cus our attention on T = properties. We have an- 
alyzed several types of DOS shapes and found three 
typical classes of the ^-dependence of the Kondo res- 
onance splitting. Since our method enables us to per- 



form NRG calculations for arbitrary band-shapes, we de- 
cide to choose an example which turns out to encom- 
pass all three classes, namely p a (oo) = ^^pD _3 / 2 (l + 
aQ)Vu~TD~T~aA, where uo <E [-D - aA,D - a A], 
Do = D + A, [a = 1(— 1) for f (j)], a square-root shape 
DOS equivalent to a parabolic band (as for free elec- 
trons) with Stoner splitting A |l4||, an d some additional 
spin asymmetry Q, which modifies the amplitude of the 
DOS [see Fig. [[{insets)]. 

In Fig. [I] we present the weighted spectral function 
A(uo) = 7re 2 /h^2 cr r cr (0)A a (uj), normalized such that for 
uo = it corresponds to the linear conductance G = A(0), 
as a function of energy uo and eo. We focus on a nar- 
row energy window around the Fermi surface where the 
Kondo resonance appears; charge resonances are visible 
when eo or U + eo approach the Fermi surface, namely 
at energies e /U > -0.1 or < -0.9. Although the NRG 
method is designed to calculate equilibrium transport, 
one can still roughly deduce, from the spin-splitting of 
the Kondo resonance of the equilibrium spectral function 
A((j), the splitting of the zero-bias anomaly AV in the 
non-equilibrium conductance G(V), since eAV ~ 2Ae 7] 
(Ae = e| — e i is the splitting of the renormalized levels). 

We present A(uo) for three DOS shapes depicted in the 
insets of Fig. d (i) Q = (a,d), (ii) Q = 0.1 (b,e), 
and (hi) Q = 0.3 (c,f), with 2 A = 0.3D Q, leading to 
the three typical behaviors. Here the parameter Q tunes 
the spin and p-h asymmetry [see the definition of p a (uj)] 
resulting in different behaviours (for a detailed discus- 
sion see the last section). For (i) we find nearly any 
eo-dependence of the spin-splitting; for (ii), a strong eo- 
dependence without compensation of the spin-splitting 
(i.e. no crossing), and for (hi) a strong eo-dependence 
with a compensation (i.e. a crossing) and a change of 
the direction of the QD magnetization. The compen- 
sation (crossing) corresponds to the very peculiar situa- 
tion where the Kondo effect (strong coupling fixed point) 
can be recovered in the presence of ferromagnetic leads 
without any external magnetic field. A behavior as pre- 
sented in Fig. da,b) was recently observed experimen- 
tally PilQji; where indeed a variation of the gate voltage 
results in two split conduction lines G(V,V g ) which are 
parallel for one case and converging for the other case, 
similar to our findings. 

Effect of a magnetic field. - In Fig. ^d,e,f) we 
show how a magnetic field B modifies the results of 
Fig.E{a,b,c): in (i) the spin splitting can be compensated 
at a particular magnetic field 5 comp (here B comp /U = 
0.017) and the Kondo effect is visible in a wide range 
of eo; for (ii), at B/U = 0.0083, the Kondo effect is re- 
covered only at one particular eo- value, which depends 
on the applied magnetic field; case (iii) shows that the 
crossing point shifts with B. Since B comp can be viewed 
as a measure of the zero-field splitting, Ae(B = 0, eo) — 
— £? comp (eo), the eo-dependence of Ae can be measured 
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FIG. 2: (color online). The QD's linear conductance G as 
a function of gate voltage eo and external magnetic field B 
for the DOS shapes (a), (b), and (c) as for Fig. da), (b), 
and (c) respectively, (d) Spin-dependent occupancy n a of 
the dot level as a function of gate voltage eo for the DOS 
shape as in Fig. d c ) an d B = 0. (d) The eo-dependence of 
the total occupancy of the dot n and magnetization m for 
the situation from Fig. d c )- ( e ) The conductance G for the 
situations from Fig. d c ~ dashed), (d - solid), and (f - long 
dashed). Parameters U, T, and T as in Fig. Q] 



by studying that of £? comp , for which one needs to mea- 
sure the linear conductance G(eo,B) as a function of 
both B and eo- In Fig. Efa-c) we plot G(eo,B) for the 
three bands of Fig. ^ The two horizontal ridges (res- 
onances) in Fig. (2fa-c) correspond to quantum charge 
fluctuations (broadened QD level) of width ~ T. The 
lines with finite slope in Fig. E{a-c) reflect the restored 
Kondo resonance and hence map out the eo-dependence 
of ^comp(eo) = — Ae(eo) when the magnetic field compen- 
sates the spin-splitting. Interestingly the spin-splitting 
and the corresponding 5 comp tend to diverge (| Ae| — ► oo) 
when approaching the charging resonance, as is best vis- 
ible in Fig. Etc). 

Such a finite slope in G(eo,B) was observed for a 
singlet-triplet transition Kondo effect in a two level QD 
(Fig. 2(d) Ref. 18]). The corresponding transition leads 
to a characteristic maximum in the valley between two 
charging resonances (Fig. 3(c) Ref. |l3|), similarly as in 
our Fig.E^e). In that system the effective spin asymme- 
try (assumed by our model) is realized by the asymmetry 
in the coupling of two QD levels []j| . 

In Fig. Etd) we show how the occupation n a and the 



magnetic moment (spin) of the QD m = n^ —n^ = 2(S Z ) 
change as a function of eo for the situation of Fig. [He). 
One finds that even though B = 0, it is possible to con- 
trol the level spin-splitting of the QD, i.e. its spin, and 
thereby change the average spin direction of the QD from 
the parallel to anti-parallel alignment w.r.t. the lead's 
magnetization. This opens the possibility of controlling 
the QD's spin state by means of a gate voltage without 
further need of an external magnetic field, which is diffi- 
cult to apply locally in practical devices. 

Perturbative analysis. - One can understand the be- 
havior presented in Fig. EJa-c) by using Haldane's scal- 
ing method J2fl| , where charge fluctuations are integrated 
out. This leads to a spin-dependent renormalization of 
the QD's level position e a and a level broadening T a . In 
contrast to Ref. [jj we consider here the case of finite 
Coulomb interactions U < oo, which means that also the 
doubly occupied state |2) is of importance. The spin- 
splitting is then given by Ae = Se^ — Sei + £?, where 
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The first term in the curly brackets corresponds to 
electron-like processes, namely charge fluctuations be- 
tween a single occupied state \o) and the empty |0) one, 
and the second term to hole-like processes, namely charge 
fluctuations between the states \a) and |2). The ampli- 
tude of the charge fluctuations is proportional to T, which 
for T >T determines the width of QD's levels. Eq. Q 
shows that Ae depends on the shape of T a (uj) for all cj, 
not only on its value at the Fermi surface. The dashed 
lines in Fig. CJa-c) show ±Ae as a function of eo [from 
Eq. (PJ] for the same set of parameters as in the NRG 
calculation, and are in good agreement with the posi- 
tion of the (split) Kondo resonances observed in the lat- 
ter. Eq. shows that the dramatic changes observed 
in Fig. Q upon changing Q are due to the modification of 
the p-h and spin asymmetry. 

Eq. (J3J predicts that even for systems with spin- 
asymmetric bands T^(co) ^ Ti(uS), the integral can give 
Ae = 0, which corresponds to a situation where the 
renormalization of e a due to electron-like processes are 
compensated by hole-like processes. An example is a sys- 
tem consisting of p-h symmetric bands, T a (uj) = T a (— lj), 
for which there is no splitting of the Kondo resonance 
(Ae = 0) for the symmetric point, eo = —U/2. For 
real systems p-h symmetric bands cannot be assumed, 
however the compensation Ae = is still possible, as in 
Fig. EJc). Eq. (J3J) also shows that the characteristic en- 
ergy scale of the spin-splitting is given by T rather than 
by the Stoner splitting A (A ^> T), since the states far 
from the Fermi surface enter Eq. (JHJ) only with a logarith- 
mic weight. However, the Stoner splitting introduces a 
strong p-h asymmetry, so it can influence the character 
of gate voltage dependence significantly. 



For a flat band T a (uj) = r a , Eq. (|3J) can be inte- 
grated analytically. For Dq ^> U, \sq\ one finds: Ae ~ 
(P T/ir) Re[ 0(e o ) -<KU + e ) ], where P = (I> -Tj/r, 
4>{x) = \I/(| + *27Fr)> an d ^( x ) denotes the digamma 
function. For T = 0, the spin-splitting is given by 



Ac-(Pr/7r)ln(|e |/|tf + eo|), 



(4) 



showing a logarithmic divergence for eo — > or U + eo — ► 
0. Since any sufficiently smooth DOS can be linearized 
around the Fermi surface, this logarithmic divergence oc- 
curs quite universally, as can be observed in log-linear 
versions (not shown) of Fig. Efa-c). For finite tempera- 
ture (T > 0) the logarithmic divergence for eo — > or 
e -> -U is cut off, Ae ~ -±Pr[#(§) + ln^], which 
is also important for temperatures T <C Tk- 

In conclusion, we demonstrated, using the extended 
NRG technique for general band shapes, the possibility 
of controlling the local exchange field and thereby the 
spin-splitting in a QD attached to ferromagnetic leads by 
means of the gate voltage. A new type of the logarithmic 
divergence of the QD's level spin-splitting was found, and 
attributed to spin-dependent charge fluctuations. 
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